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Abstract 

The negative parity bands in 220'pj^ ^^^^ analyzed within the dinuclear system model which was 
previously used for describing the alternating-parity bands in deformed actinides. The model is 
based on the assumption that cluster type shapes are produced by the motion of nuclear system in 
the mass-asymmetry coordinate. To describe the reflection-asymmetric collective modes character- 
ized by nonzero values of the intrinsic excitations of clusters are taken into account. The angular 
momentum dependence of the parity splitting and the staggering behavior of the B[E1) / B{E2) 
ratios as functions of angular momentum are explained. 

PACS numbers: 21.60.Ev,21.60.Gx 
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I. INTRODUCTION 



In the even-even isotopes of actinides and also in the heavy Ba and Ce isotopes the 
low-lying negative parity states are observed together with the usually presented collective 
positive-parity states combined into rotational or quasirotational ground-state bands. The 
formation of the positive-parity rotational or quasirotational bands are connected in general 
to the quadrupole collective motion, while the lowering of the negative-parity states is a 
signature of the presence of the reflection asymmetric collective mode. There are several 
approaches to treat the collective motion related to the reflection asymmetric degrees of 
freedom. One of them is based on the concept of the nuclear mean field [l| which has a 
static mirror asymmetric deformation or is characterized by a large amplitude of reflection 
asymmetric vibrations around the equilibrium shape. Another approach is based on the 
assumption that the reflection asymmetric shape is a consequence of the a-clustering in 
nuclei It is also known from the Nilsson-Strutinsky type calculations for light nuclei 
that nuclear configurations corresponding to the minima of the potential energy contain 
particular symmetries which are related to certain cluster structures 0, 4|. Several calcula- 
tions performed for heavy nuclei 5|-[3] have shown that configurations with large equilibrium 
quadrupole deformations and low-lying collective negative parity states are strongly related 
to clustering. We mention also a different approach to description of the properties of the 
alternating parity bands which is based on the idea of the aligned octupole phonons 8, ol 

The main idea of the cluster model developed in is that a dynamics of a 



reflection asymmetric collective motion can be treated as a collective motion of nucleons 
between two clusters or as a motion in a mass-asymmetry coordinate. Such collective motion 
simultaneously creates deformations with even and odd-multipolarities. Among different 
cluster configurations only a-cluster system — ).(^~^) (Z — 2)-|-^IIe gives a significant 
contribution to the formation of the low-lying nuclear states. Within this approach the 
existing experimental data on the angular momentum dependence of the parity splitting 
and multipole transition moments (E1,E2,E3) of the low-lying alternating parity states in 
odd and even actinides 220-2285^^^^ 223,225,227^0, 222-224,226,228-232 ^j^^ 23ipg^^ 232-234,236,238u and 
240.242pQ and the medium mass nuclei ^^^'^^^'^^^Ba, isi-issp^, i^e.usQg 153,155^^ ^^^^ i46,i48Nd 
are well described. The perfect agreement between the calculated results and experimental 
data supplies the proof of the cluster features of reflection-asymmetric states. However, 
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in our previous publications we have considered in the even-even nuclei only the low-lying 
collective negative parity states with K=0. However, there are experimental data which 
indicate a presence of the low-lying collective states related to the reflection-asymmetric 
modes which are characterized by nonzero values of K. It can be also that K is not a 
good quantum number if nuclei are located in a transitional region between deformed and 

n 

spherical ones. A good example is ^^°Th [12] whose energy spectra is a challenge for the 
theoretical approaches. 

To describe the properties of the low-lying collective states related to the reflection- 
asymmetric collective mode and characterized by nonzero values of K in the framework 
of the cluster approach we should take into account intrinsic excitations of the clusters 
forming a nucleus under consideration. The aim of the present paper is a formulation of a 
simple model which gives a quantitative description of the low-lying positive and negative 
parity bands in the spectra of nuclei like ^^'^Th. Our model is based on the assumption 
that the collective oscillations of the nuclear shape lead to the formation of the cluster- 
.ype conflgurations. This model is further development of the previously used approaches 
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The choice of the collective coordinates and the procedure of the calculation of the poten- 
tial energy and of the inertia coefficients for the Hamiltonian of the model are based on the 
concept of the dinuclear system (DNS). The concept of the dinuclear system was first intro- 
duced to explain the experimental data on deep inelastic and fusion reactions. Later on it was 
applied to the description of the nuclear structure phenomena, like alternating parity bands, 



mentioned above, and superdeformed states 
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14| . The dinuclear system {A,Z) consists 



of two fragments {Ai,Zi) and {A2,Z2) with A = Ai + A2 and Z = Zi + Z2 sticked closely 
together by a molecular-type nucleus-nucleus potential. The degrees of freedom needed to 
describe the collective excitations of such a system are related to the rotation of the DNS 
as a whole, to the relative motion of the DNS fragments, to the intrinsic excitations of the 
fragments, and to the transfer of nucleons between the DNS fragments. The latter process 
is described usually by means of the mass-asymmetry coordinate 77 = {Ai — A2)/{Ai + A2) 
or by charge-asymmetry coordinate riz = {Zi — Z2)/ {Zi + Z2). The idea that the motion in 
mass-asymmetry is responsible for the formation of refiection-asymmetric deformations of 
the medium mass and heavy nuclei is based on the observation which manifest the alpha- 
cluster dinuclear system for these nuclei. Our calculations have shown that the potential 
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energy for the alpha-cluster DNS are close or even lower than the binding energy for these 
nuclei [ll|. This is also in agreement with the fact that these nuclei are good alpha-emitters. 

Previously, we considered the mass-asymmetry motion variable as the only intrinsic col- 
lective coordinate describing DNS. However, the nuclei demonstrate yet another sign of 
reflection-asymmetry, namely, low-lying (started below 1 MeV) negative parity rotational 
bands with K 7^0 which can be related to the non axially-symmetric reflection-asymmetric 
mode. The natural way to obtain such an excitation in the framework of the dinuclear 
system model is to consider the relative angular motion of the fragments in the dinuclear 
systems illustrated in Fig. 1. It is the aim of the present investigation to extend the DNS 
model to take into account such excitations. 

Basically, there are two physically different cases which can be distinguished. If the 
heavy fragment of DNS, whose lighter fragment is an a-particle is strongly deformed, the 
light fragment is tend to stay near the pole of the heavy cluster performing the small angular 
oscillations around this position. This case which holds for the heavy actinides is considered 
in the work []. If the heavy fragment is nearly spherical, the interaction between the frag- 
ments is small to prevent the light cluster from rolling along the surface of the heavy one. 
Such a case is realized in the light isotopes of actinides, including the considered here. 



The structure of ^^°Th has been analyzed in several theoretical publications |12l . llSl . Il6 |. 
Here, we apply an approach of the dinuclear system model. 



II. MODEL 



A. Hamiltonian 



The degrees of freedom chosen to describe a system with nearly spherical heavy cluster 
are related to the rotation of the DNS as a whole, to the quadrupole oscillations of the heavy 
fragment, and to the transfer of nucleons between the fragments. The Hamiltonian of the 
model can be presented in the form 

H = Ho + V,nt, (1) 

where Hq describes independent fragments of the system and Vint describes the interaction 
between the fragments. 
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FIG. 1. Schematic picture illustrates degrees of freedom used in the model to describe the dinu- 
clear system. The orientation of the vector of the relative distance R are defined by the angles 

il{6,(p) with respect to the laboratory frame system. 

We assume that the heavy cluster is spherical and perform harmonic quadrupole oscil- 
lations around the spherically symmetric shape with frequency fiwo, while the light cluster 
stays in its ground state. Then for we have the following expression 



(2) 



where /i is the reduced mass of the DNS, R is the distance between the centers of mass of 
the fragments, fi is the operator of the number of the quadrupole phonons of heavy cluster 
and is the operator of the square of angular momentum of relative rotations of the two 
fragments 

I d . „d I 



d_ . d_ 



(3) 



Angles Q = {9,4)) (see Fig.l) describes the orientation of the relative distance vector R 
with respect to the laboratory system. 

Taking into account that for the low angular momenta the deformation of heavy fragment 
is small, the interaction between the vibrational and rotational degrees of freedom can be 
taken in the lowest order with respect to the operators of quadrupole deformation ^2/* 



d'^ + d 



(4) 



where we express the operator of quadrupole deformation through the creation and anihilla- 
tion operators of bosons 

1321, = Poidt^ + d2f.), (5) 

with /3q = ^Jh/2BujQ. 

In this work we does not consider the transfer of nucleons in the DNS exphcitely. Instead, 
we assume that the system stays in its lowest state in mass-asymmetry coordinate. The 
reason for such an assumption is that the energy of a-particle DNS is around 1 MeV higher 
than the binding energy of ^^°Th. The first excited state in mass-asymmetry will appear 
high enough to neglect its influence on the low-energy part of the spectra. Instead of 
considering the mass-asymmetry motion we take the effective value of mass-asymmetry 
r]Q in our calculations. This value, which must be at least in the interval 7]^ < tjq < 1, 
where rj^ = {A — 4) /A is the mass-asymmetry for DNS with alpha-particle, is fixed by the 
description of the lowest 1~ state. 

If we neglect the interaction term Vint in (P) the eigenfunctions of the Hamiltonian can 
be constructed as 

*(n™./0/2 = [l^rnM X Yj,]^jM) , (6) 

where \nTnj\Ii) represents the n-boson wave function of the heavy nucleus, with the seniority 
T and angular momentum Ji. Since the quadrupole oscillations have positive parity, the 
parity of the states (|6]) is determined by the angular momentum of the relative rotation of 
the fragments p = (—1)^^. The energies of states (jH]) are given in this approximation by the 
expresssion 
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nil hi 
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(7) 



The set of wave functions ([6]) can be used as a basis to construct the eigenfunction of the 
Hamiltonian H in the form 

^7M,p = J2Y1 ^nfii^hh [I'^^i^^i^l) X ^/2](/M) ' (8) 

where coefficients ai^'^'* j , should be obtained by a diagonalization of H. The matrix 
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elements of Vint between states (jH]) have the following form 



X I ' ^ I [n'r'I 'i\\{d~^ + d)\\nTli), 
h h 2 



(9) 



where the reduced matrix elements of the boson operators can be calculated using the boson 
fractional parentage coefficients 

1 1 



Vri V2/ + 1 

y/n y/2I + 1 



(10) 



B. Two-level solution 

Our numerical calculations have shown that with a good accuracy the ground-state band 
and the first excited negative parity band can be presented with a good accuracy as a 
superposition of two basis states of the form ([6]). For the ground state band the two level 
approximation yields the wave function of the form (/ can have only even values) 



=sin[7o(/)] 



I--/) X Yo 
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where 



sin[7o(/)] 
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and for the energy we obtain 
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In the last expressions A = — hcu. 

For the first excited negative parity bands we have 
.1-11-1 



, (11) 



(IM) 



vl/^-^- =sin[7i(/)] 
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where 



/ 



1/2 



sin[7i(J)] = — 



1 + 



, h ^ ( ./12" (/+l)(2J+3) 

\l ^ \\l J (21+1) I 



with the energy is given as 



1 - 



1 + 



12K)/3o\ (/ + l)(2/ + 3) 
(2/+1) 




Angular momentum / can take only odd values. 



C. Multipole Moments 



Electric multipole operators are given by the expression 



(15) 



(16) 



For the dinuclear system we assume that 

p(r) = pi(r) +p2(r), 



(17) 



(18) 



where pi (i = 1,2) are the densities of the DNS fragments. Using (ITS!) we can rewrite the 
expression for the electric multipole moments for the DNS in the following form 



Q 



A/1 



47r(2A + 1^ 



Ai,Ai+A2=A 



(2Ai + 1)!(2A2 + 1)! 



A/1 



(19) 



where 



.(A1A2) 

VAi 



-1 



1^2 



(20) 



In the last expression Q*-*^ (« = 1,2) are the intrinsic multipole moments of the DNS frag- 
ments. 

Since we assume that the light fragment is spherical and can not be excited in the con- 
sidered energy range the only nonzero moment for the first fragment is Qq^^ = Zi/v^lvr. 
The second fragment is assumed to perform the quadrupole oscillations around the spheri- 
cal shape. Thus, in the linear approximation with respect to the deformation, we have two 
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nonzero moments for the second fragment Qq = Z2I and Q2 = ^^^^^ 1^2^ ■ Thus, we 
can write the expicit expressions for the dipole and quadrupole moment of the DNS in the 
form 

Qi. = e^^i^^^i? ■ ra^) (21) 

for the dipole moment and 

Q2, = e^i^^±^i?2 . Y2,m + (22) 
for the quadrupole moment. 



D. Reduced transition probabilities 



Using expression (|8]) for the wave function and (fT9!) for the multipole operators we can 
calculate the reduced transition probabilities as the 

\<If\\Qx\\I^>\' 



B{EX,Ii I 



(23) 



2/. + 1 

The reduced matrix elements for the multipole operator Q\ between the initial state i 
and the final state j has the following form 

< IjPj\\Qx\\IiPi > 



' h h I, ^ 

I[ I2 h 

Ai A2 A 



(2A + 1) 



0A201 



(2Ai + l)!(2A2 + l)! 



^/{2i + l)(2j + 1)(2A + l)(2i2 + 1)(2A2 + 1) < 



(24) 



where Ai = A — A2, and {^}({j}) stands for the set of quantum numbers of the initial 
(final) states. 

Using the two-level solutions for the ground-state band and the first excited negative par- 
ity bands we can easily calculate the explicit expression for the intraband B[E2) transitions 
and for the B{E1) transitions between these bands. 
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In the case of the quadrupole transitions we have for the transitions between the states 
of the ground state band 

i-ql^''^ sin [7o(/)] cos [7o(/ - 2)] + gf'^y|sin [7o(/)] sin [7o(/ - 2)] 



^^(0,2) ^ (i_2) ^^^^^^j ^^^^^ _ (25) 

and for the transition between the negative parity states 
B{E2, r-P- ^ (/ - 2)"-?'-) = 



V " + ^5°'^) y^sin [7o(/)] sin [7o(/ - 2)] 

(0,2) / (2/-3)(2/ + 3)(/ + l) 
+^2 y — 2(2/- 1)(2/ + 1) — cos[7o(/)]cos[7o(/-2)]) . (26) 

In the last two expressions we have 

(2,0) AIZ2 + AjZi 2 
?o = ee// ^5 -rt , 

?2°''^ = ee//|^^ii??/3o. (27) 
For the dipole transitions between the two bands our calculations yields 
B{E1,P'- (7- l)"-f-) 



^ 27>I I ^/27^«i^ [^o(^)] [^^(^ - + V ^^^T^ [^o(^)] sin [71(1 - 1)] | , 
S(£;i, r-P- ^ (7 - 1)^--) = go' sin' [7i(^)] sin' [7o(/ - 1)], (28) 
where 

A,Z2 - A2Z, 
Qo - (^eff 



III. RESULTS OF CALCULATIONS 

The results of calculations of the energy spectra of the ground state band and of the 
two lowest negative parity bands for the '^'^^Th are presented in Fig. 2 together with the 
available experimental data. The calculation is done by the exact numerical diagonalization 
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FIG. 2. Calculated and experimental level scheme of ^^'^ Th. Experimental level, spin and parity 
assignments are taken from 

of the Hamiltonian ([1]). One can see the overall good agreement between the calculated and 
experimental spectra. The comparison of the experimental and calculated results for the 
states of the ground state and of the first excited negative parity bands are presented in 
Fig. 3 in more details. It is seen that the agreement is rather good at all values of angular 
momenta. Both, the ground-state and the first negative parity bands exhibit approximately 
an equidistant behavior which is related to the harmonic quadrupole oscillations of the heavy 
fragment. 

The calculation shows that mainly two eigenvectors of Hq are present in the wave function 
of the state of the ground state band, namely, [|||(/)) x Yq] (j^j) and [|^^^^(/ — 2)) x y2\^j^y 
The contribution of the first of them is predominant at low angular momenta. As a conse- 
quence at low angular momenta the ground state band has an equidistant spectrum with 
the energy differences determined mainly by the frequency of the harmonic quadrupole 
oscillations of the heavy fragment. With increase of the angular momentum, the distance 
between the levels is slightly increased, due to the growing admixture of the component 
[|^^^^(/ — 2)) X ^2](j^^) to the wave function. This introduces a small nonlinear depen- 
dence of the 7-transition energies as a function of the angular momentum. 
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FIG. 3. Calculated (lines) and experimental (solid circles connected by lines) energies of qround- 



state band (a) and first negative parity band (b). Experimental values are taken from 



The same equidistant structure with the frequency shghtly growing with angular mo- 
mentum holds for the first negative parity band. Again, the calculation shows that mainly 
two eigenstates of Hq are present in the wave function. Namely, [["^"^(-^ — 1)) x 
and [|^^^^(/ + 1)) X ^ij^y^j^.^j • The contribution of the later component while being small 
at /=0 is growing with angular momentum. The energy differences between the states in 
the negative parity band at low angular momentum are again determined mainly by the 
frequency of the quadrupole oscillations of the heavy fragment. 

The angular momentum dependence of the vibrational frequency uj^i, = defined as 

a half of the energy difference between the energies of two neighborhood levels of the ground 
state band and of the first excited negative parity band is illustrated in Fig. 4. One can 
see the sharp decrease of the experimental values of for the transition 10+ — t- 8+ in the 
ground state band and for the transition 13^ — ?■ 11^ in the negative parity band which can 
be a consequence of the backbending phenomena. The backbending in these bands can be 
probably related to a rotational alignment of the nucleonic orbitals as it is mentioned in 12 |. 
The model considered above does not provide a mechanism which could be responsible for 
the experimentally observed behavior of the 7-transition energies. However, it is seen from 



12 



LU 




10 12 14 16 18 



I I 

FIG. 4. Calculated (line) and experimental (solid circles connected by lines) energies of j tran- 
sitions between subsequent level of the ground state band (a) and the first negative parity band 
(b). Experimental values are taken from [12]. 

Fig. 4 that the interval of change of Umf, with angular momentum observed experimentally 
and obtained in calculations is not large. The energy differences between the neighborhood 
states in bands are varying from 185 keV up to 240 keV. By taking into account the length 
of the spectra in both energy and angular momentum, the frequency can be with good 
accuracy treated constant. 

The dependence of experimental and calculated values of parity splitting in the ground 
and the first negative parity bands, treated as a unified alternating parity band, on angular 
momentum is illustrated in Fig.5. The parity splitting is defined by the expression 12| 



s{n = E{i-) 



(29) 



2/+ 1 

It is seen from the figure that for the low angular momenta the parity splitting is positive 
and becomes negative with angular momentum increase. The possibility for the negative 
values of the parity splitting is related to the fact that the ground-state band and the first 
negative parity band are of the vibrational type. In this case the sign of the parity splitting 
is determined by the difference in energies characterizing the quadrupole vibrations of the 
heavy fragment and the rotation of the light fragment around the heavy one. Indeed, in 
the zero approximation E{I^) = — 1) + 2^^, E{I~^) = \ujl and therefore S{I^) = 
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10 12 
I 

FIG. 5. Calculated (lines) and experimental (solid squares connected by lines) values of parity 
splitting (see Eg. 129\) ). Experimental values are taken from 1121. 
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2ixRF ~ t'^ (2/+1) • '^^^S' f*^^ 2^r2 < can take negative values. In the case of a 

rotational bands, the value S{I) must stay positive, achieving a zero value for the ideal 
unperturbed rotational bands of nucleus with stable octupole deformation. The important 
feature of the spectra is an appearance at low energy of the second excited negative parity 
band which contains the states of even and odd angular momenta. The second excited 
negative parity band has an interesting features. The state with angular momentum 1=2 
is lower than the state with 1=1. With increase of angular momentum the normal level 
sequence is restored. The reason for such a behavior is related to the significant contribution 
of the states [|^^(J + 1)) x Yi]^^^^^^^ and [|^^(/+ 1)) x (j+2,m) ™® 
functions for even and odd angular momenta, respectively. In the limit if Vint is equal to 
zero, these two states are degenerate. 

In Fig.6, the values of the ratio B{E2,I — )■ / — 2)/B{E2,2^ — )■ 0^) are presented as a 
function of the initial angular momentum. As it should be in the case of harmonic quadrupole 
oscillations of the heavy fragment the values of B{E2) increase linearly with /. They does 
not show any changes in the behavior for transitions between the members of the ground- 
state band and the first negative parity band since the underlying quadrupole constituents 
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FIG. 6. Angular momentum dependence of the ration of the reduced transition probabilities 
B{E2,I — 7- (/ — 2))/B{E2,2^ — t- O"*') for the quadrupole transitions between the subsequent 
levels in the ground state band and the first negative parity band. 

in both bands are the same. Thus, the experimentally observed staggering of BE1/BE2 
ratios can be attributed to the staggering of the B{E1) values (see Fig. 7). 

Such a staggering behavior of B{E1) can be qualitatively explained with the use of 
equation ( |28|) . We can see that the reduced transition probability B(E1) for the tran- 
sition from state / of the ground state band to the state (/-I) of the first excited neg- 
ative parity band consist of two contributions, since both dipole transitions are allowed 
(see Eqs.l imiMl) ) from the component [|||(/)) x V'o]^^^^^ to [|||(/)) x 'i^i] jy.^) and from 
[|^^(/ - 2)) X ^2]^^^^^ to [|^^(/ - 2)) X Yi]^^_^^^y In the opposite case of tran- 
sition from the states of the negative parity band to the states of the positive parity be- 
longing to the ground state band, we have only one allowed transition, namely from the 
component [|||(/)) x ^i] to the component [|||(/)) x Yq^^^ The transition from 



0^^(/ + 2))xn]^ 



to 



1-2 7-2 



(/ — 2)) X F2] (J is forbidden because the dipole 
operator does not change a number of the quadrupole phonons. 

The B{E1)/B{E2) ratios as a functions of an initial angular momentum are presented 
in Fig. 7. Calculated ratios for the odd initial angular momentum (i.e. for transitions from 
the states of the negative parity) lie systematically lower than the ratios for the even initial 
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FIG. 7. B{E1) / B{E2) ratios as a function of the initial angular momentum for transitions in 
the yrast and first negative parity hands. Experimental values (filled circles) are taken from 

angular momentum (transitions from the state of the ground state band). This is in agree- 
ment with the experimental data with the exception of two data points at 13~ and 14"'" . 
As it is mentioned in 12|] the large value of B{E1)/B{E2) ratio at 13~ can be attributed to 



the loss of E2 strength in the backbending. The rather small B[E1)/B{E2) value for the 
14+ attributed to the spread of El strength due to the presence of two 13^ final states. 



IV. CONCLUSION 

We suggest a cluster interpretation of the multiple negative parity bands in ^^°Th assum- 
ing collective oscillations of nucleus in mass-asymmetry degree of freedom. This collective 
motion leads to the admixture of the very asymmetric cluster configurations to the intrinsic 
nucleus wave function and creates deformations with even- and odd-multipolarities. To take 
care of reflection asymmetric modes with nonzero values of K, the harmonic quadrupole 
oscillations of the heavy cluster is considered. The resulting energy spectrum consists of 
the ground state band and several negative parity bands which exhibit nearly equidistant 
behavior. The angular momentum dependence of the parity splitting is described. The 
possibility for the negative values of the parity splitting is related in the model to the inter- 
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play between the quadrupole vibrations of the heavy fragment and the rotational motion of 
the light fragment. We described the staggering behavior of the B{E1) / B{E2)-TaXions as a 
function of the angular momentum. The BE{1) transitions from the state of the negative 
parity to the state of the positive parity is hindered, because these states contain different 
numbers of the quadrupole phonons. The results of calculations are in overall good agree- 
ment with experimental data. This work is further development of the previously developed 



approaches 
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